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Abstract 

Networks of coupled degrade-and-fire (DF) oscillators are simple dynamical models of as¬ 
semblies of interacting self-repressing genes. For mean-field interactions, which most mathemat¬ 
ical studies have assumed so far, every trajectory must approach a periodic orbit. Moreover, 
asymptotic cluster distributions can be computed explicitly in terms of coupling intensity, and 
a massive collection of distributions collapses when this intensity passes a threshold. Here, we 
show that most of these dynamical features persist for an arbitrary coupling topology. In par¬ 
ticular, we prove that, in any system of DF oscillators for which in and out coupling weights 
balance, trajectories with reasonable firing sequences must be asymptotically periodic, and pe¬ 
riodic orbits are uniquely determined by their firing sequence. In addition to these structural 
results, illustrative examples are presented, for which the dynamics can be entirely described. 


1 Introduction 

To predict the long-term behavior in networks of interacting units is a predominant challenge 
in nonlinear science, with applications in many disciplines, from physics to biology and to the 
social sciences, to cite a few examples m- In particular, a recurrent question is to characterize 
collective properties such as synchronisation and predictability in terms of the network topology 
and interaction strengths mE]. While this problem has received considerable attention from 
theoreticians, mathematically rigorous descriptions of (global) nonlinear behavior are scarce, and 
only address limited circumstances, e.g. weak-coupling regimes m and assemblies of pulse coupled 
oscillators with excitatory coupling 01113 eg. Hence, the theory remains largely incomplete and 
network phenomenology still lacks a comprehensive rigorous footing. 

In the last years, a model for the population dynamics of simple gene oscillators was introduced 
[8], inspired from a series of experiments on colonies of synthetic genetic circuits [6lEj, and resulting 
from the simplification of more standard delay-differential equation models m- In a few words 
(see section 2 below for more details), it consists of a collection of pulse coupled oscillators with 
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inhibitory coupling [7], and is reminiscent of the well-known integrate-and-fire model in neuroscience 
|5J; however, the phenomenologies of each are distinct. 

In the case of mean field coupling, a mathematically rigorous global description of the dynam¬ 
ics, notably its clustering and asymptotic properties, was achieved for every parameter value and 
for arbitrary numbers of oscillators (and also for the continuum approximation) [21 [S]. Analogous 
features were also described for trajectories issued from typical random initial conditions [9j. In 
addition, a recent study expanded the analysis to a more elaborate model that involves a global 
activator field in the dynamics m Motivated by including more realistic features in the math¬ 
ematical analysis, the current paper aims to extend previous (deterministic) results to arbitrary 
coupling topologies on populations of arbitrary size. 

A typical property of degrade-and-fire models is their firing process (accompanied with gene 
level resetting) that occurs when the repressor field becomes (locally) negligible and can no longer 
prevent gene expression. (Instantaneous resets are used here as a naive representation of massive 
gene expression during a tiny interval of time.) In the case of mean field coupling, after cell i has 
fired, every other cell (not simultaneously firing with i ) must fire once, before i fires again. The 
ordering in which cells fire does not change from cycle to cycle (unless cells synchronize and begin 
firing together). More importantly, this periodic and exhaustive cycle of firings imposes asymptotic 
periodicity on the trajectories themselves. 

While periodic exhaustive firing may not always hold for an arbitrary coupling topology (ex¬ 
amples will be provided below), the main result of this paper (Theorem |5.1| below) states that, 
when this is the case, the trajectory must asymptotically approach a periodic configuration, pro¬ 
vided that all cells are path-wise connected through coupling and in and out weights balance at 
every node. Even though this conclusion does not a priori cover all trajectories of degrade-and-fire 
systems, it proves that a sufficiently regular firing behavior implies a regular asymptotic behavior 
of the trajectories themselves, under a mild restriction on the coupling structure. Together with 
the analysis of orbits with symmetric components, this result paves the way to a comprehensive 
understanding of the functioning of arbitrary systems of coupled DF oscillators. 

The paper is organised as follows. The DF model of Woscillators is defined in section 2 
and global well-posedness of the dynamics is proved. In section 3, we study properties of the 
firing events, and use these features to introduce non-degenerate trajectories with exhaustive firing 
sequences; such trajectories are at the center of attention in the rest of the paper. In section 4, 
we prove that there can be at most one periodic orbit associated with each such sequence, and 
provide examples of existence and non-existence, in the case of nearest neighbor coupling. Section 
5 contains Theorem 5T and its proof, while the paper is completed, in section 6, with a study of 
the full dynamics for N = 2 and N = 3 cells (assuming some coupling symmetry in the latter case). 


2 The degrade-and-fire dynamics 

We consider the degrade-and-fire (DF) dynamics of single self-repressor genes in a colony of cells, 
driven by intercellular coupling [8]. In this context, cells are indexed by {1, • • • , N} (where N G N) 
and gene expression levels at time t G R + are represented by the vector x(t) = (xi(t))fL 1 G [0,1]^. 
Intercellular coupling of genes in this population is mitigated by a repressor field Wx = {Wxi)f =1 ^ 
defined as the action of the linear operator W on the vector x, 

N 

Wxi = ^2 w ij x ji Vi G {1, • • • , N}, 

3 = 1 
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where the symbol W = {wij)f- =1 also denotes a stochastic non-negative matrix. The dynamics 
depends as well on a threshold parameter 77 G (0,1), which is assumed to be small. Finally, we 
impose that the matrix diagonal terms satisfy wu > 77 for all i. (NB: Ref. [ 8 ] assumed mean field 
coupling, viz. WX{ — (1 — e)xi + j$Y2jLi x j f° r &•) Here, we consider for now any coupling 
satisfying the condition wa > 77 for all i G {1, • • • ,7V}, and later impose additional constraints 
when asymptotic periodicity is investigated. 

With these definitions in place, the DF time evolution of gene expression levels is given by the 
following differential equation, inspired by the delay-differential equation model in m- 


&i(t) = -Sgn (xi(t)) if Wxi(t) > rj, 

(xi(t) = xi(t- 0 ) Wx .( t)<ri C 1 ) 

\Xi(t + 0 ) = 1 W Xl {t)SV- 

In other words, the dynamics in cell i consists of two phases, depending on the repressor field 
Wxi(t). 

• When Wxi(t) > 77 , the expression level Xi(t) degrades at constant speed —1, unless it has 
reached zero (in which case, it remains at zero). In this phase, if we also have xj(t) > 0 for all 
cells j such that wij > 0 (called influencing cells), the repressor level W X{[t ) also decreases 
with speed 1 . We may eventually have W X{[t ) < 77 , depending on expression level behaviors 
in influencing cells. 


• When Wxi(t ) < 77 , a firing takes place and resets the expression level to the value 1. The 
assumption wa > 77 ensures that VFx^(t + 0 ) > 77 for the repressor field in cell i after resetting. 
Hence, after every firing, the reset genes return to the degrade phase for a positive-length 
time interval. 


Accordingly, the behavior in each cell consists of an eternal succession of degrading phases inter¬ 
rupted by instantaneous firing, unless the repressor level becomes sufficiently high to prevent any 
further firing and to maintain the gene level in a vanishing stationary state. 

Prior to investigating these behaviors in more detail, we first make sure that the dynamics is 
globally well-posed. As the next statement shows, this is granted by assuming that the evolution 
begins with a degrading phase in every cell. An element x G [0,1]^ is said to be admissible if 
Wxi > 77 for all i G {1, • • • ,7V}. (NB: any x G [ 77 ,1]^ is admissible.) 

Lemma 2.1. For any admissible x G [0,1]^, equation 0 has a unique global solution such that 
x(0) = x. 

Proof. Local existence is a direct consequence of the admissibility condition. Moreover, we have 
Xi(t) = [pci —1) + for all z, provided that t > 0 is sufficiently small. In fact, this expression holds up 
until a firing occurs. 

In addition, for every solution of 0, the function t Wxi(t) is left continuous in every cell; 
hence we must have Wxi(t ) > 77 for all (z,t) (see Lemma 3.1 in [2]). Accordingly, the first time t\x 
when a firing occurs, viz. the first firing time, is given by 

t\x — inf{s > 0 : Wxi(s) = 77 for some i G {1, • • • , TV}}. 

Clearly, we have t\x < +00 (and the infimum here is actually a minimum). 

Let the firing map F be defined on admissible vectors x G [0,1]^ by Fx = x{t\x + 0). The 
assumption wu > 77 implies that Fx is also admissible. Hence the second firing time t 2 X = t\ oFx is 
also well-defined and we have Xi(t) = ( Fx{ — t + t\x) + for t G (Fx,Fx]. By induction, one obtains 


3 


an infinite sequence {tk%}keN of firing times and a unique well-defined solution on every interval 
tk-\-\x\. 

To conclude, it remains to show that lim^^ +oc tk,x — +oo. Assume for the sake of contradiction 
that = lim^ +00 t&x < +oo. By the Pidgeonhole Principle, there exists i G {1, • • • , N} and a 
subsequence { k n } ne ^ such that Wxi(tk n x) — 77 for all n and lim n ^ +00 tk n x = t^. For this i, the 
expression Xi{tk n+1 x) = (l — tk n+1 x + t^ n x) + , together with the characterisation of the firing time 
tkn +% x , implies the estimate 

rj = Wxi(t kn+1 x) > (1 - t kn+1 x + t kn x) + w u . 

Using that lim n ^ +OC) £& n+1 x — tk n x — 0, we conclude that 77 > wu : contradicting the original 
assumption on the self-influencing weights wu- □ 


3 Properties of the firing events 


As indicated above, aymptotic behaviors of expression levels depend on the repetitive nature of 
firing events and their spatial structure with respect to the population. Space-time firing patterns 
themselves hinge on intercellular coupling topology and, possibly, also on weight intensities. 

In the case of mean-field coupling, firing patterns are strongly structured: every cell i must fire 
infinitely often, and after each firing in a single cell i, every other cell j 7 ^ i must fire either before 
or simultaneously with the next firing of i (NB: the first alternative occurs for at least one cell j, 
viz. full synchronization can never occur in this system [2]). 

In other cases of coupling, DF systems may have solutions in which some genes never fire, or 
eventually stop firing. For instance, for N — 2 and 1 — W 22 > wn(> 77 ), the expression levels given 

b30 


xi(t) = 1 - 


l-o 

tl 


tl 


and X 2 (t) = 0, Vt > 0, 


define a periodic solution of equation 0 with period t\ — 1 — (which attracts every trajectory 
with initial condition (xi,X 2 ) with X 2 < xi, see section 6.1 below for a complete analysis of N — 2 
systems). 

However, by imposing that the neighbor influence be sufficiently small (i.e. weak coupling 
regime), one can make sure that any given site fires infinitely often, for every coupling topology. 


Lemma 3.1. Assume that wa > 1 — 77 for some i E {1, • • • , iV}. Then for every solution of equation 
0 with admissible initial condition x, we have xi(t) > 0 for all t G and x^(£^ + 0) = 1 for an 
infinite sequence {ti^}ke N of reset times such that 

ti i ^ and 0 t^ k-\-i ^ , V/c G M. 

wu wu 


Proof. Using that the repressor field can never be smaller than 77 in every solution [2], together 
with the weight normalisation, we have 

77 < WXi(t) < WaXift) + 1 - wu , 

which implies Xi(t) > > 0 as desired. Moreover, this inequality imposes that the length 

of every degrading phase cannot exceed 1 — ' n ~^ Wl ^ = i.e. any two consecutive resets of X{ 
have to take place within a time interval of length □ 

1 For r G R the expression |_r*J denotes the floor function of r. In particular, |_r* — 0J = r — 1 for r G N. 
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Depending on parameters, we also suspect that two consecutive firings can happen in a cell 
without the firing of any other gene in the interim. However, similarly to mean-field coupling, 
such events never occur, provided that all weights from influencing cells are sufficiently small (and 
weight conservation is assumed). 

Lemma 3.2. Assume that W is doubly-stochastic and suppose that ma Xj^Wij < for some 
i G {1, • • • , N}. Then, after any firing in cell i alone, there must be a firing in another cell j ^ i, 
before i fires again. 

Proof. We prove that, given any vector x G [0,1]^ whose (unique) maximal coordinate is X{ , there 
must be j ^ i such that Wxj < WX{. By contradiction, assume that Wx{< Wxj for all j ^ i. By 
summing these inequalities over j and using the definition of the repressor field, one gets 

N N N 

^ik x k — EE ^jk x k E<> ^ x ik) x k'> 

k =1 k =1 jj=i k= 1 

(where the equality follows from the assumption JT Wjk = 1 for all k) which in turn yields 0 < 
ZfcLiC - Nw ik )x k . 

However, letting m = maxj^ Xj and using the assumption maxj^ Wij < we get 


E(' - Nw ik )x k 

k^i 


< ^(1 - Nw ik )m = (Nwu 
k^i 


1 )ra, 


and we must have Nwu — 1 > 0. Since m < X{, we conclude that 
1 )xi, hence the inequality 0 < J2k= i(l — ^ w ik) x k is impossible. 


- Nw ik)xk < ( Nw u - 
□ 


We aim to relate asymptotic properties of trajectories to the ordering in which cells fire. To¬ 
wards that goal, for simplicity we shall focus on trajectories for which a unique cell resets at each 
firing. Indeed, when all weights are distinct, simultaneous resets in several cells are believed to be 
non-generic (measure zero) events in the long term dynamics, except for trajectories falling into 
synchrony subspaces when the coupling possesses a symmetry na. as in the mean-field case. How¬ 
ever, uniqueness in this case can always be recovered by passing to the quotient network. Therefore, 
the single cell reset limitation is barely restrictive in the analysis of asymptotic properties. 

Given such a trajectory and k G N, let i^ G {1, • • • , N} be the reset cell at kth firing. We call 
{ik}ke N a firing sequence. In this context, the notion of a trajectory compatible with a given 
firing sequence is obvious. A firing sequence in which each segment {ik}^=eN+i ^ a 

permutation of {1, • • • , } is said to be exhaustive. A trajectory in which every expression level 

vanishing at some time is reset at the next firing above is said to be non-degenerate. 

In the case of the mean-field coupling, every firing sequence must be exhaustive and every tra¬ 
jectory must be non-degenerate (unless genes cluster into groups with identical levels, as previously 
mentioned). 

For different coupling topologies, this may not always be the case. For instance, for N = 3, 

V < 

wu = w 2 i = 3??, and ®u = w 2 2 — w 32 = W 31 = 4r?, 
the expression levels given by 


X\(t) = 1 



X2 (t) = 1 — t\ 



and x 3 (t) = 0, Vt > 0, 
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where t\ — |, define a degenerate periodic solution of equation 0 in which genes fire after their 
level has reached 0, and cell 3 never fires. More precisely, we have xfit) — 0 and x%-fit) = 1 — t\ 
when Wxi(t ) = 77 for i = 1 , 2 . 

For simplicity, in the rest of this paper, we shall only consider non-degenerate solutions with 
an (eventually) exhaustive firing sequence. 

4 Analysis of periodic orbits 

In this section, we study periodic orbits that return to their initial state after every cell has fired 
exactly once. The main statement below claims that, when the graph generated by the adjacency 
matrix W T is fully connected, non-degenerate periodic orbits are entirely determined by their firing 
sequences. 

Proposition 4.1. Assume that W is doubly-stochastic and irreducible. Then, given any periodic 
exhaustive firing sequence, either no compatible periodic non-degenerate orbit exists, or such a 
trajectory is unique. 

Proof. Let l — {ik)k= 1 be the initial (fundamental) word of the firing sequence and let ( tk)^=i be 
the corresponding firing times. Consider a non-degenerate (single cell reset) trajectory with firing 
segment l for t G [0,£/v]- Its coordinates (xi k (t]y + 0))^ =1 after the Nth firing, and hence the initial 
coordinates (xi k )fi =1 if we also assume this trajectory is periodic, must be given by 

%i k — x i k (Lv + 0) = 1 — tjsf + t/e, Vfc E {1, • • • , N}. (2) 

We aim to show that there is at most one such trajectory, i.e. that the vector (£&) is uniquely 
determined by t. Notice that monotonicity of the firing times implies X{ k < Xi k+1 , and also £q > 0, 
because the trajectory is assumed non-degenerate and xi N — Xi N {tjsf + 0) = 1. (If we had £q = 0 
then i\ would fire jointly with i n and this would contradict non-degeneracy.) 

Expression a implies that the difference xi k — t^ = 1 — Lv does not depend k. To proceed, we 
separate the case 1 — £/v > 0 , for which no expression level ever reaches 0 , from the case 1 — £/v ^ 0 
when the reset cell has vanishing expression level immediately prior to its firing. (NB: this cell 
must be the only one with vanishing level, thanks again to the non-degeneracy assumption.) 

Assume the first case 1 — tjsr > 0. From expression 0 and equation (jT|) , one obtains the 
following gene expression levels immediately before the kth firing: 


j <k, 
j > k, 


hence for the repressor field, 


N 

Wx ik (t k ) = 1 - ts'E'Wi* + Wx ik ~ x ik ■ 

j=k 


The assumption = 1 f° r a H j implies ^2 k Wxi k = X& x ik- Together with the definition 

Wxi k (t k ) = rj of the firing times, this yields the following unique characterization of the firing time 


tN- 


tN = 


N( 1 - r ]) 


E N 

k= 1 2^7 = 


N 
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In particular, the inequality 1 — £/v > 0 holds iff 


N N 

N( 1 — rf) < EE w i k ij • ( 3 ) 

k =1 j=k 

We proceed similarly in the case 1 — Uv < 0. Using the non-degeneracy assumption, the expression 
levels immediately before the kth firing are now given by 

{ 1 -tk+tj if j < k , 

0 if j = k , 

1 -t N + tj-t k if j > k , 

from which we obtain 

N 

Wx ik (t k ) = 1 - w ik i k - t N W ikij + Wx ik - x ik , 

j=k +1 


and then 


Uv 


N 0- z Ul ~ Elk 

Z^/c=l 


In particular, the existence condition 1 — tw < 0 amounts to 


N N 

k= 1 j=k 

which is complementary to the existence condition ©> in the previous case. 

With tN uniquely specified in terms of the weights wij and the firing sequence i , we have that 
the firing time equations 

Wx ik (t k ) =r], k £ {1, - ■ ■ ,N} 

can be rewritten in vector form as (Id — W)x — u , where u only depends on W{j and i and belongs 
to the space 

E= : J2 Xi = 0 

orthogonal to the eigenvector (1, • • • ,1). Writing x = x\\ + x±, where x\\ = c(l, • • • , 1) T for some 
c G M and x± G E, the firing time equations become (Id — W)x± = u. Moreover, given any x±, the 
constant c is determined by using the normalization xi N = 1, i.e. c = 1 — (x±)i N . 

Now, using the irreducibility of W, the Perron-Frobenius Theorem implies that the leading 
eigenvalue 1 of W is simple, with eigenvector (1, • • • , 1) T . Therefore, the equation (Id — W)x± = u 
has at most a single solution in E. 

At this stage, the proposition is proved. However, we complement the proof by showing that 
the equation (Id — W)x± = u has a solution in E. Notice that the stochastic matrix (W N ) T is 

scrambling, i.e. for each pair of rows of (W N ) T , there is a column on which both rows have strictly 
positive entries. Indeed, the irreducibility of W implies that, for any pair z, j, there is a word 
with L < N — 2 such that 
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Using the assumption w>a > 77 for all z, we can assume that all such chains have length N. This 
shows that any two columns of W N have at least one positive element in a coincident row. That 
( W n ) T is stochastic and scrambling implies that the restriction W N |s is a contraction for the 
£ 1 -norm on R N (Exercise 2.27 in [16]). Moreover, u G E implies that the sum ^2^=0 1 W £ u G E. 
Therefore, the series 

+00 N— 1 

Y w Nk y wt ' u e ^ 

k =0 fcO 

converges in £ 1 (R N ) and solves the equation (Id — W)x± = u. □ 

For future purposes, we denote by x L the solution of the equation (Id — W)x = u in 

U L = {x = (x ik )k=i '■ 0 < x h < x i2 < ■ ■ ■ < x lN _ 1 < x iN = 1} . 

In addition to proving uniqueness, the analysis of the firing time equations can also yield 
information about the existence of non-degenerate periodic orbits (viz. whether or not the trajectory 
issued from x L is actually periodic with firing pattern z), depending on coupling topology. For 
instance, in the case of weighted global coupling 

N 

Wxi = (1 - e)xi + e'^ajxj, Mi G {1, • • • ,iV}, 

3 = 1 

where otj > 0 and J2j a j = 1 5 a periodic orbit can exist for every exhaustive firing sequence 
l — {ik}fr =1 i and it does iff e < e L where the coupling threshold w L is known explicitly [2]. (NB: 
Non-uniform weights cq 7 ^ correspond to effective repressor fields associated with clustered 
configurations of the mean field coupling.) 

For more general coupling topologies, existence may depend on the firing sequence (in addition 
to coupling intensity), as illustrated with the following example. 

Example 4.2. (Nearest neighbor couplings on periodic chains.) Let W = (wij)ij be given by 

vo 

Wij = (1 - w)Sij + — (y_ i)j + %+i)j) , Vi, j e { 1 , • • • ,N} , 

where the index 0 is identified with N, and N + 1 is identified with 1. 

(a) In the case where N = 2n is even, no non-degenerate periodic exists which is compatible with 
any firing sequence {ik}^ 7 ! 1 for which 

( ik is odd for k G {1, • • • , n} , 

[ ik is even for k G {n + 1, • • • , 2 n} . 


(b) The non-degenerate periodic orbit compatible with the firing sequence {k}^ =1 (i.e. ik — k) 
exists iff w < T 


Proof (a) From the proof of Proposition 4.1, one has that the vector u in the component equation 
(/ — W)x± — u is given by 

N 

u h = 1 - V - tN Y 1 
j=k 


w i k ij 




in the case 1 — £/v > 0, and 


N 

u ik = 1 ~ r l- w ik i k - t N w ikij 

j=k +1 


in the complementary case. For such coupling, the assumption on the firing sequence implies that, 
in both cases, the components Ui k do not depend on k G {1, • • • n} (resp. on k G {n + 1, • • • 2 n}); 
that is, we have 2 = Ui for all i. Hence the vector x±, and thus the periodic configuration x L 
itself, must satisfy the same property, viz. (x L )i +2 = (x L )i. As a consequence, all cells with even 
(resp. odd) indexes must fire together. This property violates the single cell reset assumption. 

(b) Together with the choice of the firing sequence, that the repressor field commutes with cyclic 
permutations of coordinates implies that we must have X{[t\ + 0 ) = £z+imodiv( 0 ) for all i\ this 
implies 

{x L )i = 1 - (N - i)ti, Vi € {1, • • • ,N} 

and tjsf = Nt\. Using again the symmetry of the weights W, it suffices to check existence on the 
time interval [ 0 ,U]- 

Strict monotonicity of the gene levels implies strict monotonicity of the repressor levels. Hence, 
the firing time equation Wxi(t\) = 77 implies Wxi(t) > 77 for all i > 1 and t G [0,U] as desired, as 
well as Wx\(t) > 77 for t G [0,£i). 

It remains to check that the gene levels are positive/zero as prescribed in each case 1 — Uv > 0 
or 1 — tjsi < 0. In the first case, which explicit calculations show occurs iff w < 2rj , we obviously 
have x\{t) > 0, and then X{{t ) > 0 for all i by monotonicity, for all t G [0,ti]. For such u;, the 
periodic orbit exists and its coordinates never reach 0 . 

In the complementary case w > 277 , we know that the firing cell fires from 0, i.e. x\(t\) = 0. 
Therefore, we must ensure that all other gene levels remain positive, which by monotonicity is 
implied by x^ (U) > 0. Explicit calculations show that this is equivalent to w < 2 ^T 2 1 ^ 7 ?- The 
conclusion now follows the uniqueness of periodic orbits as in Proposition |4.1| (and that 2r\ < 

□ 

5 Asymptotic periodicity of trajectories with exhaustive firing se¬ 
quences 

We are now in position to formulate the main result of this paper. 

Theorem 5.1. Assume that W is doubly stochastic, irreducible, and that there exists a pair i,j 
of cells such that WjiWy > 0. Then, for any trajectory t x(t) which is compatible with a given 
N-periodic exhaustive firing sequence with initial word t = we have 

lim x(t k N + 0) = x L , 

k —^+00 

where t£ (£ G N) is the Uh firing time of the trajectory (and as before, x L G U L is the initial condition 
of the periodic orbit associated with 1 .) 

In addition, by combining the proof below with continuity arguments, one can show that when 
the periodic orbit associated with 1 exists, for every initial condition in U L sufficiently close to x L , 
the firing sequence of the subsequent trajectory is TV-periodic with initial word this implies that 
the periodic orbit x L is locally asymptotically stable. 


9 




We note here that the assumption wywji > 0 ensures that the statement of Theorem |5.1| holds 
for any exhaustive firing sequence t , regardless of the coupling strength (and in particular, whether 
or not sites fire after reaching zero along the orbit x(t)). This assumption can be relaxed if one 


after the proof. 


assumes instead a weak coupling regime, i.e. wu > 1 — 77 for all i G {!,••• , A^}; see Remark 5.3 


Proof. Given x(0) = x G U L , let Fh x = x{t\ + 0 ) be the reset configuration after the first firing. 
Assuming that the trajectory x(t) is non-degenerate, we have (Fi 1 x)i > 0 for all i, and so the map 
Fi x itself is defined by 


(Tq x)i 


Xi-ti if i^iu 
1 if i = i \. 


The reset maps F{ k for k G {2, • • • , N} are defined similarly. In order to prove the Theorem, we 
consider the return map 


F n — F o 

r i ~ r lN U 


0 o Fj jl , 


and its orbits that never leave the open set U L . The map F^ is continuous and piecewise affine. 
We are going to prove that a sufficiently high iterate (F^) k is a contraction (for the i 1 norm). 
It follows that every orbit in U L mu st a pproach the unique fixed point x L of i 7 ^, which must be 
located in the closure U L . Theorem 5.1 immediately follows. (NB: The periodic orbit associated 
with l exists iff x G U b ). 

In order to prove the promised contraction, we first introduce a change of variables. Let A = 
Id — IT denote the graph Laplacian operator on M. N . The expression of the firing time £&, and hence 
the expression of the image Fi k x , depends on whether the level xi k is reset from a positive value, 
or from 0. Explicit calculations yield 


Xi 


' l k 


tk 


T] - Ax k 

y-Axk 

I-™** 


if x. 
if x, 


'I'k 


't'k 


tk> 0 , 
tk < 0 , 


recalling that since W was assumed irreducible, we have w kk < 1 for all A; G {1, • • • , iV}. This 
expression implies the following commutation relation 


A o Fi k — Gi k o A , 


where the maps Gi (i G {1, • • • , N}) are also continuous and piecewise affine, with linear parts 
depending on whether X{ is reset from a positive value (+ sign) or from 0 (— sign). The linear 
parts Li^± are given by 


Li ? -\- Xj 


Xj + WjiXi if j ± i, 
waXi if 3 = G 


and Li-Xj — 


Xj + 


w jj 

l—wa 


Xi 


0 


^ j rr 

if j = i, 


for all j G {1, • • • , N}. Note that the L^± s are column-stochastic by the double stochasticity of 
W. 

The change of variable x 1 —> Ax has the following properties: 

• A is injective on U L , for every permutation l of {1, • • • , N} (see end of proof of Proposition 
[tT] above), 

• A(R n ) C E = {x G R n : E { Xi = 0}. 
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Consequently, we have 


(F*) k = A" 1 o (Gf ) k o A, Vfc G N, 

where = Gi N o- • -oG^oGq, and we regard A -1 as a map £ —»• Z4- The contraction of (F^) k for k 
large enough follows from the following statement, together with the fact that ||A - 1 ||i||A||i < + 00 . 

Lemma 5.2. There exist k G N and 7 < 1 that 

11(77* - (Gf r )7||i < 7lk - y||i, V*, y € E. 

Proof of the Lemma. The map Gf^ is continuous and piecewise affine with 2 ^ pieces, each corre¬ 
sponding to a choice of sign {+, —} for the firing event at each site in {1, • • • , N}. Labeling each 
piece by a symbolic word s = (sk)^ =1 G {+, — } N , we write 


— Li N , SN o 


o L 


1 


Because the firing map is continuous across these piecewise domains (see Proposition 5.4 in [ 2 ]), to 
prove the statement of Lemma 5.2 it will suffice to show the existence of k for which all matrices 
(L i:S ) k are cont ractions on £ in the ^-norm. To this end (similarly to the end of the proof of 
Proposition 4.1), it is enough to show that the row-stochastic matrices ((L^ s ) fc ) T are all scrambling. 
However, for every i G {1, • • • , A}, the matrix entries ( L^±) nrn satisfy 


(f'i,—)nm ^ 0 -^ (Z/2,+)m n ^ 0 • 

Therefore, one only has to prove that ( L k ) T is scrambling for k sufficiently large, where L = L l ^_^n. 

Up to a relabeling of the cells we can assume that i appears before j in the word l , viz. we 
have i = %k x and j = i^ 2 with k\ < & 2 . We are going to prove the existence of k G N such that 


{L k ) ih > 0, Vji e {1, • • • ,N} , (4) 

which is merely a restatement of the scrambling property for (L k ) T . To show Q, we use the 
following expressions for the entries of the composing matrices 

(Li,_) ifc = — ^—S ik , VMj€{l,- ,N} and (!*,_)# = 1 - , Vj G {1, • • • ,N} . (5) 

J- Wa 

By irreducibility of W, given an arbitrary j\ G {1, • • • , A}, let (jk)k =1 (where L < N — 1) be the 
shortest word such that 

W ijL W jL3L -1 • • • W J2J1 > 0 • 

Let k\ G {1, • • • , A} be such that — ji- By © and since 1 — Wi kl i kl < 1 , we have 

‘ ’ ' I J i 1 ,-)j 2 ji P w hh • 

(This is a consequence of the relations = 1 and (Li k • • • Li %f F)j 2 j 1 = 

(Lj 1 -)j 2 j 1 .) Let now &2 be such that ik 2 = J 2 and consider separately the cases £q < k^ and 
k 2 < k\. In the first case, we have 


and then 


‘ ’ * ^4 1 + i,-)j3j2 — (I J ik 2i —)jsj 2 P. W jsj2 ? 


feor ’ ’ ’ —)j3ji — w j: 


3332 W 323i ’ 
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In the second case, we have 


0 , 

so no positive estimate holds for Lj 3 j 1 . However, we certainly have 

Lj2jl — W hh anC ^ i^k 2 - ' ' ' ^H,-)j3j2 = (^ / 4 2 ,-)j3j2 — W j3j2 > 

hence 

(^j k 2 — ’ * * — ^)j3jl — W j3j2 W j2jl • 

By repeating this process, we obtain that there exists n G {0, • • • , L — 1} for every (jk)k=i such 
that 

(Lik L - • • • I J i 1 -L n )ij 1 > ^ij L Wj L j L _ 1 • • • Wj 2 j 3 > 0 . 

To conclude, we use the following property 


C^hv>— ’ ’ ’ I J ik L + 1 ,—)« 1 if f 0 , • • • , £jv }5 

> WjiWy if t G {ife L+ i, • • • , Lv}, 

where the second estimate follows from the fact that j appears after i in l. Using w^wy > 0, it 
results that (L nJrl ) i j l > 0 for every {jk}k=v 

Finally, this last estimate also implies La > w^wy > 0. Hence, letting k = maxj -^l ^ n, we can 
always multiply by L k ~( n+1 ^ to obtain the desired estimate Q. □ 

Remark 5.3. In the weak coupling regime wa > 1 — r], the map G ^ is affine with linear part 
L l ^ + \n. Using the properties 


— ^jk( 1 &ij) T 

instead of one can repeat the proof mutatis mutantis, now using the estimate 

’ ’ ’ ^^k L + 1 j +)» — ^ 

in the case i G {ik L + 1 , • * * , £jv}; to obtain the conclusion ©• 

6 Low dimensional examples 

In complement to previous results on arbitrary DF systems, we now present examples for which 
the dynamics is entirely known. Ignoring the case of mean-field coupling which has previously been 
described for populations of arbitrary size TV, we focus on low dimensional systems TV = 2 and 
TV = 3. 

6.1 Two coupled oscillators (TV = 2) 

Letting W 2 = wu and w\ = W 21 for simplicity, we get 


Wxi = (1 - w Z -i)xi + W3-iX 3 -i, i = 1,2, 

and we consider separately the cases w\ + W 2 < 1, w\ + W 2 = 1 and > 1. 

• Case w\ + W 2 — 1. We have VFxi = Wx 2 for all x G [0, l] 2 ; hence cells 1 and 2 must fire 
simultaneously and evolve in sync after the first firing, as an TV = 1 DF oscillator. 
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Figure 1: Trajectories in the unit square [0, l ] 2 for systems of two coupled DF oscillators. Single 
arrows represent motion at speed 1. Double arrows represent resets at firings. Left. w\ + W 2 < 1: 
every firing sequence must be exhaustive with cells 1 and 2 firing alternatively. Right. W 1 +W 2 > 1: 
Depending on its initial location with respect to the diagonal x\ — X 2 , any trajectory reaches in 
finite time a periodic orbit with either x\ — 0, or X 2 — 0. NB: The white (un-shaded) region 
of [0, l ] 2 corresponds to the set of admissible initial conditions. Recall that we always assume 
wu = 1 — ws-i > rj for i = 1, 2, so that the dynamics is well-defined. Moreover, W{ < 77 corresponds 
to the weak coupling case, when expression level Xi cannot reach 0 (see Lemma 3.1). Conversely, 
when Wi > 77 , the level X{ may vanish, depending upon the trajectory. 


Of note, this property extends to any population size TV e N in the trivial case when no weight 
Wij depends on i, viz. WX{ — ^2jWjXj for all i. Similarly, full synchrony in the trajectory holds 
for any coupling W when the initial coordinates X{ do not depend on i. For convenience, from here 
on, we shall assume that not all weights, nor all coordinates, are equal. 

• Case W 1 +W 2 < 1 . Until it fires, any trajectory initially located in the segment x\ — 1 lies below 
the diagonal x\ — X 2 of the unit square [0, l] 2 , and cell 2 must fire first (see left panel in Fig. [l]). 
Conversely, in any trajectory initially located in the segment X 2 — 1, cell 1 must fire first. It results 
that every firing sequence (of a trajectory out of the diagonal) must be exhaustive, and cells 1 and 
2 must fire alternatively. 

If we also assume that w\ = W 2 > 0 (so that W is doubly stochastic and irreducible since we 
also assume 1 — > 77 from the beginning), using Theorem 5.1, we conclude that every trajectory 

asymptotically approaches the periodic orbit associated with the firing pattern { 1 , 2 }, for which 
the corresponding periodic orbit always exists. 

• Case w\ + W 2 > 1. As opposed to the previous case, for any trajectory initially located in the 
segment x\ = 1, only cell 1 can fire; hence X 2 eventually reaches and stays at 0 . All such trajectories 
reach in finite time the periodic orbit where cell 1 oscillates alone (whose expression is given before 

: right panel in Fig. [I]). Similarly, any trajectory initially at X 2 — 1 reaches a periodic 


Lemma 3.1 


see : 


trajectory with x\(t) = 0 . 


13 

























6.2 Three coupled oscillators (TV = 3) 

For simplicity, we assume that each cell influences all other cells in the same way, W 21 = w 3 i (= w±), 
w \2 = ws 2 {= ^ 2 ) and W 13 = W 2 s(= 1 C 3 ), i.e. we have 

{ VFxi = (1 - W 2 ~ W 3 )x 1 + 1C 2 X2 + W 3 X 3 
WX2 = + (1 - Wi - W$)X2 + W 3 X 3 

Wx 3 = W1X1 + W2X2 + (1 -w>2~ W 3 )x 3 

Thanks to this symmetry, for all pairs i, j of indexes, we have Wxi = Wxj in the plane X{ — x j, 
and then Wx\ — Wx 2 = Wx 3 along the diagonal of the cube [0, l] 3 . In order to characterise the 
dynamics, we need to determine which parts of the planes Wxi — 77 can be reached under the flow. 
Similarly to TV — 2, we shall separate the cases w\ + W 2 + w 3 < 1 and w\ -\- W 2 ~\~ w 3 > 1. 




Figure 2: Locations of firing planes WX{ — 77 in the unit cube [0, l ] 3 for systems of three coupled 
DF oscillators with weights satisfying W 21 — w 3 i(= irq), W 12 — w 3 2 {— ^ 2 ) and — W 2 3 {— w 3 ). 
Left. W 1 +W 2 +W 3 > 1: In the pyramidal cone with section delimited by the square x\ — 1 and with 
apex at the origin, the plane Wx\ — r\ (green shaded area) lies above the other planes WX 2 — rj 
and Wx 3 = 77. Since this cone contains all trajectories starting with x\ — 1, cell 1 must fire first, 
and by induction, this cell is the only one to fire in this case. By symmetry, similar arguments 
apply to the cones respectively delimited by X 2 = 1 and x 3 = 1. Right. w\ + W 2 + ^3 > 1: Now, 
the pyramidal cone x\ — max^rj decomposes into two regions. In the lower sector X 2 > x 3 , the 
plane Wx 3 = r] (blue shaded area) lies above the two other planes, while the plane Wx 2 — r\ (red 
shaded area) dominates in the sector x 2 < x 3 . By symmetry, we conclude that the firing sequence 
of every non-degenerate trajectory must be exhaustive and periodic with pattern { 1 , 2, 3}, {1, 3, 2}, 
or one of their cyclic permutations. 

• Case w\ + W2 + ^3 > 1 . Any trajectory initially located in the square x\ — 1 remains in the 
pyramidal cone x\ — max^xj (delimited by the planes x 3 — 0, x\ = x\ — x 3 and x 2 = 0), until 
it fires. In this sector, the plane Wx\ — rj lies above the two other planes WX2 — 77 and Wx 3 — r\ 
(more precisely, the truncated solids respectively delimited by the planes Wx2 — r\ and Wx 3 = 77 
in this sector, both contain the truncated solid delimited by the plane W x\ — rj, see Fig. [2j left 
panel), hence cell 1 fires first and is reset to the square x\ = 1. By induction, it results that the 
trajectory remains in this cone forever. Moreover, both X2 and x 3 must eventually vanish and the 
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trajectory reaches in finite time, a periodic orbit with X 2 (t) = x%(t) = 0. Similar scenarios occur 
for trajectories initially located in, respectively, the squares X 2 = 1 and x% — 1 . 

Of note, the phase portrait in this case is fully preserved under asymmetric perturbations of 
parameters, provided that the following conditions hold: 


1 - wsi — ws2 < wi 3 A W 23 , 1 W21 W23 < w\2 A W32, and 1 - w\2 — rei3 < W21 A 11731. 


• Case wi + W 2 + < 1. In this case, the pyramidal cone x\ — max^J decomposes into two 

regions, according to the sign of X 2 — £3 (ignoring the case X 2 = £ 3 ). In the lower sector X 2 > £ 3 , 
the plane Wx% — 77 lies above the two other ones (Fig. [2j right panel), while in the upper sector 
%2 < ^ 3 , the plane Wx 2 — r\ dominates. 

Moreover, a trajectory initially in the lower sector is located, after the first firing that resets 
cell 3, in the sector x\ > X 2 of the square £3 = 1. Therefore, cell 2 must fire second, then cell 1, 
and then the trajectory is back into the lower section of x\ — 1. By induction, it results that the 
trajectory has exhaustive firing sequence with periodic pattern {3,2, 1 }. As known from Theorem 
5.1, in the doubly stochastic and irreducible case w\ — W 2 — W 3 > 0 (which is equivalent to 


mean-field coupling), it must asymptotically approach the associated periodic trajectory. 

Similarly, trajectories starting in the upper sector X 2 < x% of the pyramidal cone x\ — max^{a;J 
have exhaustive firing sequence with periodic pattern {2,3,1}. Moreover, the fate of trajectories 
starting form the other cones X 2 — max^xj and X 3 — max^J can be obtained in the same way, 
by applying the permutation symmetries. Hence, the dynamics is also fully described in this case. 
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